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The multiconfiguration Dirac-Hartree-Fock (MCDHF) method was employed to calculate atomic 
electric dipole moments (EDM) of the superheavy element copernicium (Cn, Z = 112). The EDM 
enhancement factors of Cn, here calculated for the first time, are about one order of magnitude 
larger than those of Hg. The exponential dependence of the enhancement factors on the atomic 
number Z along group 12 of the periodic table was derived from the EDMs of the entire homologous 
series, ioZn, ^igCd, ^ioHg, ifiCn, and iHUhb. These results show that superheavy elements with 
sufficiently long half-lives are potential candidates for EDM searches. 
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I. INTRODUCTION 


The existence of a non-zero permanent electric dipole 
moment (EDM) of an elementary particle, or in a nonde¬ 
generate system of particles, would be one manifestation 
of violation of parity (P) and time reversal (T) symme¬ 
tries Qi- Violation of P symmetry has been observed in 
the /3-decay of ®°Go i followed by decay of muons Q and 
pions Violation of charge and parity (CP) symme¬ 
try has been observed in the weak decay of neutral kaons 
. Violation of T symmetry is in turn equivalent to 
violation of combined CP symmetry, through the com¬ 
bined CPT symmetry, which is considered invariant Q. 
Both CP and T symmetry violations have been observed 
in the neutral kaon system Q, although direct T sym¬ 
metry violation has been disputed [1, [13] ■ More recently 
a direct observation of the T symmetry violation in the 
B meson system has been reported [Il|. The violations 
of P, C, CP, and T symmetries are predicted by the 
Standard Model (SM) of particle physics [3 in. How¬ 
ever, the SM leaves several issues unexplained, such as 
the origin of baryogenesis, the mass hierarchy of funda¬ 
mental particles, the number of particle generations, the 
matter-antimatter-asymmetry observed in the universe, 
and the nature of the dark matter. These and other is¬ 
sues are addressed within a large number of extensions of 
the present version of the SM. Several of these extensions 
predict EDMs induced by the P and T violating inter¬ 
actions and also EDMs of the fundamental particles sig¬ 
nificantly larger than the values predicted by the SM it¬ 
self. The predictions can be tested, and searches for per¬ 
manent electric dipole moment are underway presently 
in various systems — neutrons [l^ . electrons in para- 
and diamagnetic atoms [3) [I3| , molecules [Hill , and 
other species The experimental searches have 


not yet detected a non-zero EDM, but they continue to 
improve the limits on EDMs of individual elementary 
particles, as well as limits on CP-violating interactions, 
parametrized by the coupling constants Ct and Cp (see 
references Hfj, for details, and the Table II in the 
reference [23| for a summary). 

The primary objective of the present paper is the cal¬ 
culation of EDM for the superheavy element coperni¬ 
cium We evaluated the contributions to the 

atomic EDM induced by four mechanisms [l3|- tensor- 
pseudotensor (TPT) and pseudoscalar-scalar (PSS) in¬ 
teractions, nuclear Schiff moment (NSM), and electron 
electric dipole moment interaction with the nuclear mag¬ 
netic field (eEDM). In each case we show that there is 
an order of magnitude increase of atomic EDM between 
mercury and copernicium. The second objective of the 
present paper is to derive the Z-dependence of atomic 
EDM. We show that numerical EDM results are consis¬ 
tent with an exponential Z-dependence along the group 
12 elements. 


II. MCDHF THEORY 

We used the MGDHF approach to generate numeri¬ 
cal representations of atomic wave functions. An atomic 
state function Ali^PJMj) is obtained as a linear combi¬ 
nation of configuration state functions ^{jrPJMj) that 
are eigenfunctions of the parity P, and total angular mo¬ 
mentum operators and Mj: 

^'(yPJMj) = ^ cMurPJMj). (1) 

r 

The multiconfiguration energy functional was based on 
the Dirac-Goulomb Hamiltonian, given (in a.u.) by 
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where a and /3 are the Dirac matrices, and p is the mo¬ 
mentum operator. The electrostatic electron-nucleus in¬ 
teraction, V(rj), was generated from a nuclear charge 
density distribution p (r), which was approximate d by 
the normalized to Z two-component Fermi function |26l| 

= i + e(%/a ’ ( 3 ) 

where a and b depend on the mass of the isotope. The 
effects of the Breit interaction, as well as QED effects, 
were neglected, since they are expected to be small at the 
level of accuracy attainable in the present calculations. 


III. MCDHF WAVE FUNCTIONS 

We calculated the wave functions of five diamag¬ 
netic atoms of group 12, and subsequently the EDMs 
in the ground states of the entire homologous series, 
3 oZn, “|Cd, ^soHg, ?i 2 Cn, and fgiUhb. The numer¬ 
ical representations of the wave functions were gen¬ 
erated with the relativistic atomic structure package 
GRASP2K [ 23 , based on the multico nfig uration Dirac- 
Hartree-Fock (MCDHF) method [ 2 ^ 128143^ . Electron 
correlation effects were evaluated with methods described 
in our previous papers (ssl - l^ . Core-valence and valence- 
valence correlations were included by allowing single and 
restricted double substitutions to five sets of virtual or¬ 
bitals. The full description of numerical methods, vir¬ 
tual orbital sets, electron substitutions, and other details 
of the computations can be found in [s^. Compared 
with [s^ the double electron substitutions were however 
extended from the nsnp to the (n— l)dnsnp shells in the 
present paper (see section below for details). 


IV. EDM CALCULATIONS 


An atomic EDM can be written as a sum over states 
(equation (4) in reference [1^): 




En — Ei 


( 4 ) 


where | 0 ) represents the ground state | 4 '( 7 PJMj)) of a 
closed-shell atom from the group 12, with J = 0 and 
even parity. The summation runs over excited states 
|4'(7i(—P)JiAfjJ), with Ji = I and odd parity. A cal¬ 
culation of an atomic EDM requires evaluation of the 
matrix elements of the static dipole D^., and the matrix 
elements of the Hint interactions, which induce an EDM 
in an atom [s^- In order to perform these calculations 
the CRASP2K package was extended. The extension in¬ 
cludes programs for matrix element calculations, based 
on spin-angular integration [s^- Here Hint represents 
one of the four interactions mentioned in the section |l] 
above, Eq and Ei are energies of ground and excited 


states, respectively. The full description of the EDM the¬ 
ory underlying the presents calculations can be found in 
reference [33|, and in references therein. 

The summation in equation o involves an infinite 
number of bound states, as well as contributions from 
the continuum spectrum. The sum over the bound spec¬ 
trum was evaluated by explicitly calculating contribu¬ 
tions from the lowest five odd states of each symmetry us¬ 
ing numerical wave functions. Then the method of ’Rie- 
mann zeta tail’, described in reference [s^ . was applied 
to sum up the contribution from the remaining bound 
states. To this end we showed that a summation over a 
Rydberg series, when extrapolated to large values of the 
principal quantum number n of the running electron (and 
where the energy denominator saturates at the ionisation 
energy) converges to the Riemann zeta function. The ex¬ 
plicit numerical summation accounts for 98 percent of the 
whole sum, and we evaluated the upper bound on the rest 
(the infinite tail) of the sum by exploiting regularities of 
the Rydberg series. The relative correction, i.e. the to¬ 
tal contribution from the trailing terms (called Riemann 
zeta tail) divided by the total contribution from the five 
leading terms, is of the order of 1.5 percent for mercury 
i99Hg, below 2 percent for copernicium ^®®Cn. We 
neglected the Riemann zeta tail correction for the other 
three elements (®®Zn, ^^^Cd, ^®^Uhb). 

The contribution from the continuum is difficult to es- 
tim ate, since it is partially accounted for by the virtual 
set [38j In the present paper we computed only the con¬ 
tribution of the bound states. We neglected the explicit 
summation over continuum, and assumed that the con¬ 
tinuum spectrum contribution were included into the er¬ 
ror budget. The evaluation of the sum over the contin¬ 
uum part of the spectrum could in principle be carried 
out either fully numerically [s^, or again by an extrap¬ 
olation, based on the fact, that the oscillator strength 
density is continuous across the ionization threshold , 
and above mentioned regularities carry over to the con¬ 
tinuum spectrum. 

The electronic matrix elements in equation (|4]) are not 
isotope-specific. However, the atomic wave functions do 
exhibit a (rather weak) dependence on the atomic mass of 
the element of interest, through the nuclear electrostatic 
potential, which depends on the nuclear charge density 
distribution, which in turn depends on the nuclear mass 
number, through the equation ([3]). Therefore, all numer¬ 
ical results in Tables DiniiiHi and m were obtained for 
specific isotopes, such as ^®®Hg and ^®®Cn, and they do 
exhibit a (negligibly weak) dependence on atomic masses. 


V. MERCURY 

The calculations for ^®®Hg were performed in a simi¬ 
lar manner as those presented in our previous paper . 
The results from DF and from calculations with the first 
two layers of virtual orbitals (i.e. the first three lines in 
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TABLE I: TPT interaction contributions to EDM in different 
virtual sets, in units (ca) |e|cm), for ®®Zn, 

and ^®®Cn, compared with data from other methods. 
See text for explanations and details. 



'5S>Zn 

i“Cd 

199 Hg 


VOS Th SE 

Th SE 

Th SE 

Th Th2 Th3 


DF -0.07 -0.07 -0.35 -0.36 -7.29 -6.15 -59.86 -61.50 -66.66 

1 -0.08 -0.09 -0.39 -0.45 -4.13 -4.86 -48.53 -50.95 -53.95 

2 -0.09 -0.11 -0.45 -0.54 -4.66 -5.23 -58.38 -58.92 -62.96 

3 -0.10 -0.12 -0.47 -0.57 -4.84 -5.53 -59.31 -64.53 -68.76 

4 -0.10 -0.12 -0.48 -0.59 -4.79 -5.64 -57.67 -61.04 -65.26 

5 -0.11 -0.12 -0.49 -0.60 -4.84 -5.64 -57.51 -60.75 -64.98 


Ref. f4lUDHFl 

-2.4 

Ref. f42UDHFl 

-2.0 

Ref. f4lUCI4-MBPTl 

-5.12 

Ref. f41im,PA'l 

-5.89 

Ref. f42URPA'l 

-6.0 

Ref. [43UCPHF1 

-6.75 

Ref. [44](CCSD) 

-4.3 


Tables mini mi an d IIVI) are in fact identical with the re¬ 
sults published in |^. Further calculations differ in the 
scope of the double electron substitutions, which were 
extended from 6s6p to 5d6s6p shells. The results of the 

TABLE II: PSS interaction contributions to EDM in different 
virtual sets, in units (l0~^^Cp {ua) |e|cm), for ®®Zn, ^^^Cd, 
and ®®®Cn, compared with data from other methods. 
See text for explanations and details. 



“iCd 

199Hg 



VOS Th SE Th SE 

Th SE 

Th Th2 

Th3 


DF-0.13-0.14-0.94-0.96-25.47-21.49-199.52-252.66-274.11 

1 -0.15-0.17-1.05-1.21 -14.54-17.16-199.52-209.13-221.73 

2 -0.19-0.23-1.19-1.46-16.38-18.39-240.22-242.15-259.07 

3 -0.20-0.24-1.25-1.53-17.01-19.47-244.96-266.95-284.65 

4 -0.20-0.24-1.28-1.58-16.84-19.84-237.56-251.33-268.95 

5 -0.22-0.24-1.30-1.60-17.02-19.85 -236.88-250.07-267.78 

Ref. [^(DHF) -8.7 

Ref. [^(CI+MBPT) -18.4 

Ref. [41]('RPAi _ -20.7 _ 

calculations are presented in Tables HI HIllIIIl and IIVI The 
number of virtual orbital sets (VOS) is listed in the first 
column of each table (see chapter III of reference [s^ for 
definitions and for the details of the calculations). The 
line marked ’DF’ (Dirac-Fock) in the VOS column rep¬ 
resents the lowest-order approximation, with zero sets 
of virtual orbitals. It should be noted that the values 
in the Tables marked ’DF’ and ’DHF’ are not equiva¬ 
lent. Those marked ’DF’ were obtained in the present 
calculations with only the two lowest excited states in¬ 
cluded in the summation in equation (|T|). The results 
marked ’DHF’, obtained with MBPT methods, involved 
summation over the entire spectrum of virtual orbitals, 
using various methods to construct the virtual orbital 
set ^,112,113 • Neither ’DF’ nor ’DHF’ include electron 
correlation effects and therefore they are relevant only 


TABLE III: Schiff moment contributions to atomic EDM in 
different virtual sets, in units {l0“^^[S/(|e| fm®)] |e|cm}, for 
®®Zn, ^^^Cd, and ®®®Cn, compared with data from 

other methods. See text for explanations and details. 




69- 

Zn 

iiiCd 

199 

Hg 


98'! Cn 


vos' 

Th 

SE 

Th SE 

Th 

SE 

Th 

Th2 

Th3 

DF 

-0.04 

-0.04 

-0.18 -0.19 

-2.86 

-2.46 

-17.73 

-17.26 

-19.53 

1 

-0.05 

-0.06 

-0.21 -0.26 

-1.95 

-2.45 

-13.64 

-12.96 

-14.53 

2 

-0.06 

-0.07 

-0.25 -0.32 

-2.11 

-2.42 

-17.05 

-15.96 

-17.78 

3 

-0.06 

-0.08 

-0.27 -0.34 

-2.21 

-2.58 

-20.09 

-22.66 

-24.58 

4 

-0.06 

-0.08 

-0.28 -0.35 

-2.19 

-2.62 

-17.75 

-18.02 

-19.95 

5 

-0.07 

-0.08 

-0.28 -0.35 

-2.22 

-2.63 

-17.62 

-17.77 

-19.71 


Ref. f4lUDHFl 

-1.2 

Ref. f4lUCI+MBPTl 

-2.63 

Ref. f4lURPA') 

-2.99 

Ref. [^(CI+MBPT) 

-2.8 

Ref. f46UTDHF) 

-2.97 

Ref. f44UCCSD~) 

-5.07 


TABLE IV: Contributions of electron EDM interaction with 
magnetic field of nucleus, to atomic EDM in different virtual 
sets, in units {de x 10“^), for ®®Zn, ^^^Cd, ^®®Hg, and ®®®Cn, 
compared with data from other methods. See text for expla¬ 
nations and details. 



"’’Zn 

“’Cd 

199Hg 


VOS Th SE 

Th SE 

Th SE 

Th Th2 Th3 


DF 

0.13 

0.14 

-0.62 

-0.63 

16.04 

13.41 

314.03 

324.40 

350.09 

1 

0.11 

0.09 

-0.64 

-0.71 

8.47 

9.58 

254.78 

269.22 

283.51 

2 

0.13 

0.13 

-0.69 

-0.81 

9.63 

10.64 

305.55 

309.48 

328.86 

3 

0.14 

0.14 

-0.72 ■ 

-0.85 

9.99 

11.30 

305.13 

329.18 

349.47 

4 

0.14 

0.14 

-0.73 

-0.87 

9.90 

11.53 

300.39 

318.41 

338.60 

5 

0.13 

0.11 

-0.75 ■ 

-0.88 

10.00 

11.50 

299.67 

317.11 

337.40 


Ref. [411(DHF) 

4.9 

Ref. [£3 (DHF) 

5.1 

Ref. [4£(CI+MBPT) 

10.7 

Ref. [4t](RPA) 

12.3 

Ref. [471fRPA) 
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for the purpose of evaluating the contributions of elec¬ 
tron correlation for the expectation values of interest. A 
larger number of VOS represents in principle a better ap¬ 
proximations of the wave function. The line marked ’5’ in 
the VOS column represents the final approximation, with 
five sets of virtual orbitals (MCDHF-VOS.5, represented 
by red circles in Figure[2|). The difference between VOS.4 
and VOS.5 may (cautiously) be taken as an indication of 
accuracy. For each element the calculated values of the 
energy denominators in equation (|T]) were used to evalu¬ 
ate the atomic EDMs. These fully theoretical EDM val¬ 
ues are marked ’Th’ in Tables HI HIl Hill and IIVI Semiem- 
pirical EDM values (marked ’SE’ in the Tables) were also 
evaluated for ®®Zn, ^^^Cd, and ^®®Hg, with the energy 
denominators taken from the NIST database [d^l. Level 
identifications were made with the atomic state functions 
transformed from jj-coupling to LSJ coupling scheme, 
using the methods developed in [H, • 
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VI. COPERNICIUM 

Three different sets of energy denominators for ^®®Cn 
were used. Those from the present calculations are 
marked ’Th’. For comparison purposes we computed also 
the EDMs with the energy denominators taken from two 
other theoretical papers [5l|, . The results in column 

marked ’Th2’ were obtained with the energy denomina¬ 
tors taken from [5l|, who used a large-scale MCDHF 
method. The authors of [5l| evaluated also the ioniza¬ 
tion limit of copernicium and their calculated ionization 
energy was used in our evaluation of EDMs for those 
levels which were not reported in M- The energy de¬ 
nominators in the columns marked ’Th3’ were taken from 
reference , where the energy spectrum was computed 
with two methods: CI-I-MBPT and CL We gave priority 
to the CI-I-MBPT results; the Cl results were used when 
CI-I-MBPT data were not available; for the remaining 
levels the energy denominators were replaced by the cal¬ 
culated ground state ionization energy. The accuracy of 
our calculated energy values, as well as those from the 
references [5l| and , is better than 20 % for the lowest 
excited levels of mercury. 

The mass number 285 for the element Cn was chosen 
due to predictions that heavier isotopes are more sta¬ 
ble than the lighter ones . The lifetimes of sev¬ 

eral known isotopes of Cn are counted in minutes [s^, 
which make them amenable to atom traps, and subse¬ 
quent spectroscopy. It is predicted that still heavier iso¬ 
topes of Cn, with mass numbers in the range 290-294, 
may have half-lives counted in years [H^. 

We observed a similar pattern of contributions from in¬ 
dividual electronic states, as described in [^. The triplet 
6snp and the singlet 6snp ^Pi states are the dominant 
contributors to atomic EDM in the ^®®Hg spectrum. For 
the ^®®Cn case the dominant contributions arise from the 
lowest states of symmetries, i.e. Isnp^Pi, 7snp^Pi, 
Altogether they contribute in excess of 98% of the total 
EDM. The remaining Rydberg states contribute less than 
2 percent. Instead of an explicit error analysis for the cal¬ 
culations of EDM for ^®®Cn we applied a comparison with 
mercury. Estimates of the magnitudes of EDMs induced 
by the TPT, PSS, NSM, and eEDM mechanisms in mer- 
cury, have been performed with several theoretical meth¬ 
ods [ill, lil, 1^ . With one or two exceptions [3, [sll , 
they all agree within reasonable error bounds — of the 
order of 10-20 percent [3|. The results of the MCDHF 
calculations for mercury, both in present paper as well as 
in [3]: are well within these bounds. We expect that the 
present calculations for ^®®Cn, performed with the same 
MCDHF model as those for ^®®Hg, would also fit within 
error bounds of similar size. 


VII. UNHEXBIUM 

In addition to the calculations described above we 
have done uncorrelated DF calculations for fg^Uhb and 


for iBe. There are several theoretical predictions [ 13 - 
[3| which suggest that the heaviest homologue in the 
Zn-Cd-Hg-Cn-Uxx group would not be element E162 
(Unhexbium), but E164 (Unhexquadium). Due to a 
very large spin-orbit splitting of the 8p shell, the rel¬ 
ativistic 8 pi /2 shell becomes occupied before the 7d 
shell is filled [ 13 ]. Therefore, at the end of transi¬ 
tion metals in the row eight of the periodic table ap¬ 
pears the element E164, with the ground configura¬ 
tion [Cn]5g^^6f^^7d^^7p^8s^8p^, with all inner shells 
closed, and with two electrons in the 8 pi /2 shell (the 
8 pi /2 shell is, in fact, also closed). However, the pres¬ 
ence of the 8p shell would complicate the calculations 
of EDMs, and, more importantly, it would complicate 
comparisons along the homologous series. Therefore 
we have deliberately chosen a (doubly artificial) isotope 
fg|Uhb, of element E162, with electron configuration 
[Cn]55i®6/i'‘7di°7/8s2. 


VIII. Z-DEPENDENCE 


Atomic properties depend in various ways on the 
atomic number Z, both in isoelectronic sequences [60- 
I3| . as well as along homologous sequences [ 3 , 1 ^ . 
In many cases approximate analytic relations were de¬ 
rived [ 3 13 - 13] , and several atomic observables exhibit 
a polynomial or power dependencies on the atomic num¬ 
ber Z. 

Atomic enhancement factors of the PT-odd interac¬ 
tions in neutral atoms scale with nuclear charge as dat ^ 
c?Z^. The factor Z^ arises from an estimate of the 
strength of the electric field in the vicinity of an atomic 
nucleus (see chapter 6.2 of the reference [l|), but it has 
been pointed out that on top of the Z^ enhancement 
of the PT-odd interaction there is another enhancement 
factor, arising from relativistic contraction of the elec¬ 
tronic wave function P, [3j7lj] . 


Kr 


r(3) y 

T{2j + 1)) [ ao J 


( 5 ) 


Z-dependence of atomic EDMs induced by the (P, T)- 
odd Hint interactions is governed by the Z-dependence 
of three factors in equation (|T|): matrix element of the 
(P, r)-odd Hint operator, matrix element of the elec¬ 
tric dipole Dz operator, and the energy denominator 
{Eq — Ei). The matrix elements of the electric dipole 
operator are constrained by the Thomas-Reiche-Kuhn 
rule. In case of the elements of group 12 they are fur¬ 
ther constrained by the Wigner-Kirkwood sum rule (see 
chapter 14 of the reference [3])- The two lines, ns^ ^S'o — 
nsnp ®Pi and ns^ — nsnp ^Pi, dominate the Wigner- 
Kirkwood sum in all five elements, making the matrix 
element of approximately constant along the homolo¬ 
gous series. Transition energy denominators in the equa¬ 
tion ffl do not depend on Z along the homologous se¬ 
ries [ 3 ]! except for small variations due to shell con- 
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tractions, shell rearrangements, etc (excluding the Uhb 
element, with its large spin-orbit splitting mentioned in 
the section Ivll abovel. 

Therefore, the dominant role in establishing the Z- 
dependence of atomic EDMs along the homologous se¬ 
quence is taken by the operators. Following the 

analysis in chapter 8 of the reference [l[ , it can be shown 
that in the vicinity of a point-like atomic nucleus the large 
PriK and small Quk radial components of a one-electron 
wave function may be expressed as 


Pn>i{r) 


QnK^^) 



7) 




(^ 

2 

/ ao y-'>' 


) r (27 + i) 

V2Zr/ 

(6) 


2 / 

' ao y-y 


r (27 + i) 1 

.2Zr) 


(7) 

where k is the angular momentum quantum number, 
7^ = — a^Z'^, a is the fine structure constant, is the 

effective principal quantum number, and oq is the Bohr 
radius. The radial integrals involved in the calculations 
of matrix elements of Hint include the integrands of the 
combinations of the large and small radial com¬ 
ponents, of the type {PaPb ± QaQb) or {PaQb ± PbQa)- 
All these integrals include factors in the integrands which 
effectively cut off the integrals outside the nucleus [s^, 
and eventually Z-dependence of the atomic EDM in the 
form 


dat 


Vr(27+ 1 )/ V ao j 


( 8 ) 


is obtained, where k depends on a particular form of 
the integrand and where r^v is the effective cut off ra¬ 
dius. The right hand side of the equation ([8]) is dis¬ 
played in Fig. [TJ All four combinations {PaPb, QaQb, 
PaQb, and PbQa) of the one-electron wave function fac¬ 
tors from equations ([6]) and © are represented as func¬ 
tions of atomic number Z. The index a represents nsi /2 
orbitals, the index b represents np ^/2 orbitals. The nu¬ 
clear radius has been computed using rjq = 
where rg = 1.25 fm. The relation of atomic mass A to 
atomic number Z has been evaluated from the neutrons 
to protons ratio NjZ = 1-1- l‘2aA, derived from 

the Bethe-Weizsacker formula [t^, with ac = 0.711 and 
a A = 23.7. The empty circles in the Fig.[T]show positions 
of the four elements (Zn, Cd, Hg, Cn) considered in this 
paper. Neglecting a weak Z-dependence through the 
gamma function 2/r(27-|-1), for small values of Z the 
polynomial factor Z^ determines the functional form of 
the dependence on Z while for large values of Z the ex¬ 
ponential factor (2ZrAr)^^“^ takes over. It can be shown 
numerically, as can bee seen in the Figure [U that the 
polynomial Z^ shape dominates up to about Z = 60, 
then in the range 60 < Z < 120 the function dat{Z) is 
approximately exponential, and eventually the approxi¬ 
mation breaks down, because the analytic approximation 



FIG. 1: (Color online): Z-dependence of the atomic EDM. 
The right hand side of the equation ([Sll , calculated from (abso¬ 
lute values of) one-electron wave function factors PaPb, QaQb, 
PaQb, and PbQa- The factors were generated from the equa¬ 
tions and 0, and evaluated at r = rjv, as functions of 
atomic number Z. See text for details. 


in equations [Bland [7] is valid only within the atomic num¬ 
ber Z range, where bound solutions of the Dirac equation 
exist (Z < 137 for point-like nuclei). 


The analysis above has been made under the assump¬ 
tion of a point-like Coulomb field in the Dirac radial equa¬ 
tion. The finite sizes of nuclei entered only when the 
equation (l8|) was evaluated. For extended nuclei the so¬ 
lution of the Dirac equation depends on the specific form 
of the nuclear charge distribution IlMlTi. Flambaum 
and Ginges 0, and Dzuba et al |2l| assumed uniform 
distribution of the electric char ge i nside a sphere (with 
the normalization factors from [TOj), and obtained en¬ 
hancement factors of a similar form as in equations (|5|) 
and (O, for angular symmetries Si/ 2 , P 1 / 2 , and P3/2. 


In the present paper the numerical calculations for 
the homologous series of the group 12 of the periodic 
table (sgZn, ^^gCd, ^IgHg, n|Cn, and igiUhb) were 
performed with extended nuclear model m, for which 
bound solutions of the Dirac-Fock equations exist up to 
Z = 173 [7^. The dependence of EDMs on atomic num¬ 
ber Z along group 12 of the periodic table is presented 
in Fig. [T] The red circles represent our final results, cal¬ 
culated within the MCDHF-VOS.5 electron correlation 
model described above. The blue pluses represent the 
uncorrelated DF results. The green plus in the upper 
right corner represents the EDM value for Uhb. Due to 
very large spin-orbit splitting of the 8p shell, (see sec¬ 
tion |Vl] above), the Uhb energy denominators are dis¬ 
tinctively different from those of other members of the 
homologous series. To compensate for this splitting, we 
also computed the EDMs for Uhb with energy denomi¬ 
nators taken from Cn. The latter value is represented by 
the square in Fig. [T] The solid line is fitted to the four 
(Zn, Cd, Hg, Cn) final results. The dashed line is fit- 
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FIG. 2: (Color online): Atomic EDM (absolute values) in¬ 
duced by the NSM as a function of atomic number Z. Red 
circles = MCDHF-VOS.5 results with 5 virtual orbital sets. 
Blue pluses = uncorrelated DF results (0 sets). The lines are 
exponential functions fitted to the four points, representing 
Zn, Cd, Hg, and Cn. Solid red line = fit to MCDHF-VOS.5 
results. Dashed blue line = fit to uncorrelated DF results. 
The lines are extrapolated beyond Z = 112. The two sym¬ 
bols in the upper right corner represent Uhb (excluded from 
the fitting). Green plus = DF result for Uhb with calcu¬ 
lated Uhb energy denominators. Green square = DF result 
for Uhb with energy denominators taken from Cn. The sizes 
of circles represent approximately the relative accuracy of the 
MCDHF-VOS.5 calculations. See text for details. 


ted to the four uncorrelated DF results. The Uhb results 
were excluded from the fitting. The regression analysis 
yields the following relations: 


^TPT 

[- 1 . 22 ( 8 )-e°°™®(®)'^ 

- 5(6)] • 10-22 

d^^^ = 

[ - 30(l)-e° °8i3(3)-z 

-8.54(1)] • 10-26 

^NSM _ 

[—1.77(7) 

-k 2 ( 2 )] • 10-^6 

d^EDM = 

[ 2.74(8)-e°-°®'^^(2)'^ 

- 15(9)] • 10-6 


where the numbers in parentheses represent RMSE devi¬ 
ations. The third line of the equation © is displayed in 
Fig.d 

Similar regression analysis can be done for the semi- 
analytic relations presented for the point-nucleus case in 
the equation (|S]) and in the Fig. [H but restricted to the 
range of atomic numbers 30 < Z < 112, covered by the 
four elements (Zn, Cd, Hg, Cn) considered in this paper. 
The analysis yields ^ g 0 .oi 7 -z ^o. 022 -z^ 

somewhat smaller exponents than those presented in the 
equation ([HI). 

The deviation of the EDM value for the element E162 
from the fitted function in Figure [2] may be explained 
by several possible mechanisms: rearrangements of the 
valence shells, i.e. relativistic contraction of the 8 s and 
8 pi /2 shells, which results in the above mentioned large 
spin-orbit splitting of the 8p shell; variation of transition 
energy denominators, induced by shell rearrangements; 
contribution of QED effects, which could be quite size¬ 
able near the end of the periodic table at Z = 173 


However, the most likely explanation is the breakdown 
of the exponential approximation near the end of the 
periodic table. The analytic approximation in equa¬ 
tions [S] and [7] is valid only within the atomic number 
Z range, where bound solutions of the Dirac equation 
exist ( Z < 137 for point-like nuclei, Z < 173 for ex¬ 
tended nuclei). The element E162 is close to the end of 
the periodic table at Z = 173, where determination of 
a numerical wave function, even at the Dirac-Fock level, 
may be problematic or impossible, and one might ex¬ 
pect a question whether perturbation theory still works 
in QED for elements close to Z = 173 [t^. 

At very short distances Z-dependence algebra is dom¬ 
inated by the cutoff radii rjv (related to the sizes of the 
nuclei), and by the power series solutions for PnK and 
Quk. at the origin [32l |. The power series coefficients for 
PriK and Qnn depend on the nuclear potential (again re¬ 
lated to the sizes of the nuclei), and are constrained by 
orthogonality of the one-electron functions with the same 
symmetry. The dominant contributions to the matrix el¬ 
ements of the Hint operators come from the valence ns^ 
orbitals in the ground state, and from the lowest npi /2 
and np 3/2 orbitals in the excited states. 

The upper graph in the Fig. [3] shows the coefficient 
Po of the lowest order polynomial in the series expan¬ 
sion at the origin of the large component P of the radial 
function of the valence orbitals (ns, npi/ 2 , np 3 / 2 ) of the 
elements from the group 12 (plus beryllium). The quan¬ 
tum number n assumes the values 2, 4, 5, 6 , 7, 8 for 
Be, Zn, Cd, Hg, Cn, Uhb, respectively. The lower graph 
in the Fig. [3] shows the atomic EDMs induced by the 
TFT, PSS, NSM, and eEDM mechanisms, as functions 
of atomic number Z for the elements of the group 12 
(plus Be). For the purpose of this comparison, all values 
in Fig. [3] were obtained in the Dirac-Fock approximation, 
without account of electron correlation effects, and with 
the extended nuclear model ([3]). Analogously to the point 
nucleus case (represented by the equation ([3])), the func¬ 
tion dat{Z) is approximately exponential in the range 
60 < Z < 120, i.e. for heavy and superheavy elements 
relevant from the point of view of the EDM searches. 
Both graphs in Fig. [3] show similar Z-dependencies as 
those in Figure[Tl i. .e. the polynomial shape up to about 
Z = 60, then approximately exponential in the range 
60 < Z < 120, and eventually the exponential approx¬ 
imation breaks down near the end of the extended Pe¬ 
riodic Table of Elements [H, [t^ at Z = 173. When 
comparing the shapes of the curves in the upper and 
lower graphs, one has to bear in mind that radial inte¬ 
grals in matrix elements of the Hint operators involve 
valence ns^ orbitals in the ground state, and npi /2 and 
np ^/2 orbitals in the excited states. The apparent simi¬ 
larity of the npi /2 and np ^/2 curves in the upper graph 
and the four curves in the lower graph is a numerical 
confirmation of the dominant role of power series coef¬ 
ficients in the matrix element of the Hint operators, as 
well as of the proportionality relations between matrix el¬ 
ements, established in [2l| . Beryllium does not belong to 
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FIG. 3: (Color online) Upper: power series coefficients po 
of valence orbitals as functions of atomic number Z. Blue 
squares = np^/ 2 \ red diamonds = ns\ green circles = npi/ 2 \ 
n = 2, 4, 5, 6, 7, 8 for Be, Zn, Cd, Hg, Cn, Uhb, respectively. 
Lower: atomic EDM (in arbitrary units on logarithmic scale) 
induced by NSM (red squares), TPT (blue circles), eEDM 
(magenta stars), and by PSS (green diamonds), as a function 
of atomic number Z. All lines in both graphs are drawn only 
for the guidance of the eyes. See text for details. 


the group 12 (which results in the visible deviation of Be 
from the fitted function) but was included in Fig.[3]to in¬ 
dicate that the dominant role of power series coefficients, 
as well as the proportionality relations 21|, are not lim¬ 
ited to one group of elements. The other deviations from 
linearity in the Figure [3] are induced by electron corre¬ 
lation effects, whose contributions differ from element to 
element due to shell rearrangements. 


IX. CONCLUSIONS 

Electric dipole moments (EDMs) have not yet been de¬ 
tected experimentally. The experimental searches have 
been going on for the last 50 years, and the role of the¬ 
ory is not only to provide the limits on the fundamen¬ 
tal parameters, but also to guide the experimentalists to 
atoms, molecules, and other systems with suitable en¬ 
hancement factors. Experimentalists need to know (or¬ 
der of magnitude of) enhancement factors before they set 
up their apparatus to detect EDM in a new species 
The present paper is intended to present the calculations 
of EDMs, carried out with the multiconfiguration Dirac- 
Hartree-Fock theory, of a superheavy element 2f|Cn. The 
main conclusion of the present paper is the suggestion 
for setting up an EDM experiment on a superheavy el¬ 
ement, which would result in an order of magnitude in¬ 
crease of sensitivity, compared to a homologous heavy 
element. Such homologous pairs include (but are not 
limited to) Yb-No, Hg-Cn, T1-E113, Po-Lv, At-E117, 
Rn-E118, Fr-E119, Ra-E120. If the exponential Z- 
dependence derived in the present paper is assumed for 
all above homologous pairs, an increase of sensitivity 
by a factor 8-30 should be expected. The best limit 
on the EDM of a diamagnetic atom comes from ^®®Hg, 
for which d(^^ |^Hg ) < 3.1 x 10“^®e-cm (95% C.L.) has 
been reported Our calculations indicate that for 

the Hg-Cn pair the increase of sensitivity would be 
57.5/4.8, 236.9/17.0, 17.6/2.2, and 299.7/10.0 for TPT, 
PSS, NSM, and eEDM, respectively. Over the last 50 
years the precision of EDM experiments has been improv¬ 
ing by about an order of magnitude per decade |2l.ll8ll78l- 
l8l| . On this timescale an experiment on Cn would be 
equivalent to time travel into the future over a distance 
of about ten to twenty years. 

We are of course aware of the fact that an EDM ex¬ 
periment on a short-lived superheavy element is imprac¬ 
tical at this time. However, the techniques for trappin g 
atoms [H, , controlling quantum systems [U , 

and performing spectroscopic investigations of radioac¬ 
tive [86| and superheavy elements [87| advance rapidly. 
At the same time the quest for the superheavy island 
of stability continues [Hj, and sooner or later one 
may expect a breakthrough of laser spectroscopic meth¬ 
ods into the realm of superheavy elements [s^. 

The EDM experiments with superheavy elements, if 
ever becoming feasible, would probably constitute the fi¬ 
nal frontier for atomic tests of violation of parity (P) and 
time reversal (T) symmetries. In recents years the molec¬ 
ular avenue promises to become more competitive in the 
EDM searches. The advantages of molecular eEDM ex¬ 
periments is in the large values of the effective electric 
field, several orders of magnitude higher than those in 
atoms Current progress in cooling and trap¬ 
ping molecules Mil , as well as molecular ions [9^, , 

may soon allow to increase coherence times and im¬ 
prove population control in molecular experiments, which 
might translate into a significant advantage of molecular 

















experiments over atomic ones. 

While an EDM experiment on a short-lived superheavy 
element is impractical at this time, still less practical 
would be an experiment on a superheavy molecule. How¬ 
ever, when molecular EDM experiments become feasible, 
one may also envisage making, trapping, and eventually 
performing spectroscopy of superheavy molecules. It is 
difficult to say what is impossible, for the dream of yester¬ 
day is the hope of today and the reality of tomorrow ■ 
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